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ABSTRACT 



We present sub-arcsec angular resolution, high-Strehl ratio mid-IR adaptive 
optics images of the powerful OH/IR source and cool hypergiant NML Cyg at 
8.8, 9.8 and 11.7 fim. These images reveal once more the complexity in the 
dusty envelope surrounding this star. We spatially resolve the physical structures 
(radius ~ 0."14, ~ 240 AU adopting a distance of 1.74 kpc) responsible for NML 
Cyg's deep 10 /im silicate dust absorption feature. We also detect an asymmetric 
excess, at separations of ~0."3 to 0."5 (~520 to 870 AU), NW from the star. The 
colors of this excess are consistent with thermal emission of hot, optically thin 
dust. This excess is oriented in the direction of the Cyg 0B2 stellar association, 
and is likely due to the disruption of NML Cyg's dusty wind with the near-UV 
radiation flux from the massive hot stars withi n Cyg B2. This interaction was 
predicted in our previous paper ISchuster et al\ (l2006al ). to explain the geometry 
of an inverted photo-dissociation region observed at optical wavelengths. 



Subject headings: stars: individual(NML Cyg) — circumstellar matter — stars: 
mass loss — open clusters and associations: individual (Cyg 0B2) — stars: 
imaging — techniques: high angular resolution 
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1. Introduction 

The powerful OH/IR source and cool hypergiant NML Cyg is o ne of the most mas- 



sive and luminous M stars i n the Galaxy ("^40 Mp^, M} 



2006a 



Schuster 



2003) 



-9.5. iHv: 



Humphreys &: Davidson 1979 : Morris fc JuralE983 : Humphreys fc Davidson 



and et al 



1994 



1972 



Schuster et al 



Originally discovered by iNeugebauer. Martz &: LeightonI (119651 ). 
NML Cyg was easily identified as an extremely bright infrared source (/ ~ 8, ~ 1 
mag). The star' s relatively close proximity to the Cyg 0B2 stellar association confirms its 
high luminosity (IMorris fc Jural Il983l : ISchuster et al\ l2006al ) . Even though NML Cyg is lu- 
minous, it is also heavily obscured at visual wavelengths due to high interstellar (IS) and 
circumstellar (CS) extinction, with V fainter th an 16.6 mag. Its visual/near-IR spectrum 
indicat es a M6 spectral type (IWing et al\ 19671 ). and its spectrum peaks in the 2-20 /im 
range Jcillett et al\ Il968l : Istein et all Il969f )~ Thus, the spectrum indicates a substantial. 



optically thick dusty CS envelope obscuring the central star. NML Cyg's thick envelope is a 
result of its strong post-Main Sequence (post-MS) w ind (with a velocity of 23 km s~^) and 
high continuous mas s- loss rate of 6.4 x 10~ ^ Mq yr~^ ( Hyland et al. 1972; Bowers et a/.lll98ll : 
Morris fc Jura 1983 ). Habing et al. (1982) mapped an unusual oblong shaped HII regioE0 at 
21 cm nearly 1' in size around NML Cyg, revealing another unique feature of this enigmatic 
M supergiant. In addition, more recent high- angular resolution, high-contrast images from 
the Hubble Space Tele scope (HST) WFPC2 camera reveal an asymmetric CS nebula of dust 
scattered stellar light (ISchuster et a/.ll2006al ). The CS nebula has a shape similar to the HII 
region, but is about 300 times smaller in radius. It is necessary to consider NML Cyg's local 
environment and its effects on the star in order to understand this object's true nature. 

Although NML Cyg is a very luminous post-MS star which is losing mass at a prodigious 
rate, it does not currently dominate its local interstellar environment because of its unique 
location within the Milky Way. It lies in relatively close proximity to Cyg 0B2, which is 
possibly the largest OB stella r association - in size, mass and number of OB stars - in the 
Galaxy (see iKnodlsederl |2003| . and references therein for a multi- wavelength review). The 
Cyg 0B2 asso ciation spans nearly 2° on the sky, or ~30 pc in radius at the distance of 
1.74 ±0.2 kpc (IMassey fc Thompson! 1 199 ll ). making it one of the closest massive associations 
to the Sun. Cyg B2 has ~120 O st ar members, including 5 of the 10 most luminous O 



stars in the Galaxy (iHumphreysl Il978l ) . and approximately 2500 B stars. The stellar mass 



within the association is possibly as high as 10^ Mq. The O stars provide the UV radiation 
responsible for a Stromgren sphere that extends to a radius of at least 2.° 74 on the sky, or 
>80 pc (the projected separation between Cyg 0B2 and NML Cyg). The Lyman continuum 



^See references therein for the history of the discovery. 
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photon flux of tlie wliole association is estimated to be 10^^ pli s~^, or >10^ pli cm~^s~^ at 
tfie location of NML Cyg. Mid-IR images from the Midcourse Space eXperiment (MSX) 
satellite show NML Cy g is located in a relative void of gas/dust inside the Cygnus X super 
bubble (ISchusterl 120071 ) ■ This configuration allows the Lyman continuum photons and near- 
UV radiation from the hot stars within Cyg 0B2 to travel the large distance to NML Cyg 
relatively unimpeded. 



Morris fc Jural (119831 ) demonstrated that the Lyman continuum radiation from Cyg 0B2 



is responsible for the asymmetric HII region observed around NML Cyg. The HII region 
is the result of a photo-ionization interaction between the spherically symmetric, expanding 
hydrogen wind from NML Cyg balanced against an incident plane par allel Lyman continuum 
photon flux from Cyg 0B2 (see Figures 1 and 2 in lMorris fc Juralll983l ). The photo-ionization 
is inverted in the sense that usually the nearest massive star is the ioii ization source, and 
the ionization occurs outward from the central source. Morris fc Jura also demonstrated 



that the strength of the Lyman continuum flux from Cyg 0B2 and the density of atomic 
hydrogen around NML Cyg are sufficient to produce the observed 21 cm emission, thus 
confirming NML Cyg's high continuous mass-loss rate. In a similar way, ISchuster et al 



(l2006al ) described how near-UV photons with en ergies between a few and 13.6 eV photo- 
dissociate the molecular gas in NML Cyg's wind. ISchuster et a/.l 's model explains why the 
asymmetric shape of the CS nebula resolved in Hubble Space Telescope WFPC2 images is a 
consequence of this physical interaction with Cyg 0B2. Recent observa tions of excited-state 
OH maser emission in NML Cyg's wind by ISjouwerman et al\ (120071 ) are consistent with 
an increased OH column density, presumably arising from dissociated H2O molecules at the 
water photo-dissociation region. Since the atomic and molecular gas around NML Cyg is 
disrupted, the dust facing Cyg 0B2 is unprotected from the near-UV photons. It is thus 
likely that the CS dust grains are heated and destroyed by the radiation from Cyg 0B2 and 
should emit a detectable infrared signature. 

In Section[2|, we present sub-arcsec angular resolution, high Strehl Ratio MIRAC3/BLINC 
mid-IR images of NML Cyg obtained at the 6.5 meter MMT observatory with the Adap- 
tive Optics (AO) system. These observations are among the highest resolution images of 
this star at mid-IR wavelengths, and reveal the complexity in the asymmetric dust envelope 
surrounding NML Cyg. New photometry from the Spitzer Space Telescope's Infrared Array 
Camera, the Hubble Space Telescope's WFPC2 camera, and ground-based infrared observa- 
tions are presented in Section [2l The Spectral Energy Distribution (SED) for NML Cyg's 
CS envelope presented in Section [3] shows an optically thick silicate dust absorption feature. 
In Section HI we characterize the resolved structures in the dusty CS envelope. In Section 
we show that the CS dust component facing Cyg 0B2 exhibits a silicate feature in emission. 
We suggest that UV radiation from the hot, massive stars within Cyg 0B2 likely causes the 
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inversion in the silicate feature through external heating and dust destruction. We conclude 
by discussing these results in relation to previous works on NML Cyg, and comparing NML 
Cyg to the extremely luminous M-type hypergiant VY CMa. 



Observations 



2.1. MMT/AO Observations 

The inverse Photo-Dissocia tion Region (PDR) model for NML Cyg's CS envelope de- 
scribed in lSchuster et al\ (l2006al ) implies the presence of warm dust near the photo-dissociation 
front (s). According to this model, Cyg 0B2's UV radiation is heating/destroying the CS 
dust on the side facing Cyg 0B2. In order to corroborate this model through direct ob- 
servation, we obtained high-angular resolution, high-contrast mid-IR images of NML Cyg's 
dusty CS nebu la. The observati ons were made using the University of Arizona's adaptive 
optics system Jwildi et a/.l l2003h and the MIRAC3/BLINC camera Jfloffmann et a/.lll998 



Hinz et a/.ll2000l ) on the 6.5 meter MMT telescope. Table [T] summarizes the observations 



made in July 2006 and Figure [T] shows the NML Cyg images at 8.8, 9.8 and 11.7 /im. The 
BLINC module was used in 'chop' mode with a frequency of ~1 Hz, which is useful for back- 
ground (sky) subtraction when combined with 'nod' dithers. Using natural guide star mode, 
the MMT/AO system produced stable, high-Strehl ratio (~96 — 98%), high signal-to- noise, 
nearly diffraction limited Point-Spread Functions (PSF) at 10 /xm. We have used the source 
itself as guide star for the AO system (while extended in the infrared, NML Cyg is a point 
source at optical wavelengths where the AO system operates). With a magnitude R ~ 14, 
NML Cyg is at the sensitivity limit for the MMT/AO wavefront sensor, but was successfully 
used thanks to the excellent sky conditions. 

Figure [T] also shows images of the IR standard 7 Dra, which was observed at two 
separate epochs and sky orientations for calibration and PSF characterization (see Table 
[H only one epoch is shown in Figure [1]). The PSF stability was excellent throughout the 
night and the sky conditions were especially good, with seeing <0."25. In these PSF images 
the minima/maxima in the diffraction pattern are closely matched to the Airy pattern for a 
6.35m apertur(fl, with a FWHM of 0."33 at 8.8 /im, 0."36 at 9.8 /im, and 0."44 at 11.7 /xm. 

The nearly diffraction limited images have similar benefits in resolution and stability 
to images from space-based observatories, and at mid-IR wavelengths we can detect and 



^The MMT adaptive secondary is undersized to improve IR performance, resulting in a 6.35m equivalent 
resolution limit. 
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separate the externally heated dust from the rest of NML Cyg's CS nebula. Relatively 
long exposure times and a 5 point dither pattern, repeated multiple times with sub-pixel 
shifts, provided very high dynamic range images of NML Cyg and the calibration/PSF 
standards, in the range ~ 2 x 10^ — 10^ (PSF peak relative to background noise). The images 
were processed with the Drizzle algorithm t o maximize the quality of the data and recove r 
sampling resolution ( Hook fc Fruchteiill997l : iKoekemoer et a/.ll2000l : iFruchter fc Hookll2002l ). 
The processed MIRAC3 images were resampled by 2x2 from 0."0954 pix~^ to a scale of 
0."0477 pix-i. 

The 8.8, 9.8 and 11.7 /im filters were chosen to provide good spectral coverage of 
the 10 /xm silicate dust feature observed for this luminous cool star (IMonnier et all 119971 : 
Blocker et a/.ll200ll ). The photometry of the source in the three MIRAC3 bands is listed 
in Table [2J The flux in each band was measured in an aperture of 2" radius, with the 
sky residual background flux (after sky subtraction for each chop/nod sequence) determined 
outside this aperture. Photometric calibration was obtained using the PSF reference star 
7 Dra as photometric standard. The photometry has been corrected to a common airmass. 
The uncertainty in the photometry quoted in Table [2] do not include 5 - 10% systematic 
uncertainties typical for ground based mid-IR observations. 



2.2. Spitzer /IUAC Observations 



Table [T] summarizes our ob servations of NM L Cyg obtained in July 2004 with the 
Infrared Array Camera (IRAC, iFazio et al\ |200J) on-board the Spitzer Space Telescope 
JWerner et a/.ll2004j : ICehrz et adl2007h The IRAC images (GTO Program 124, AOR 6588416, 
pipeline version S14.0.0) were processed with IRACproc, a softw are package that facilitates 
co-addition and analysis of IRAC data (ISchuster et a/.l l2006bl ). IRACproc also improves 
the identification and removal of cosmic rays and other transient signals from the processed 
images, while simultaneously improving image quality, photometric accuracy and sensitivity 
(signal-to-noise). Due to the extreme brightness of the source, the core of the IRAC im- 
ages was saturated, preventing the adoption of standard aperture photometry techniques to 
measure the flux in the IRAC bands. By fitting a High Dynamic Range (HDR) PSF to the 
diffraction spikes and extended "wings" of the IRAC images, we have however derived the 
source photometry in the IRAC bands at 3.6, 4.5 and 8.0 /zm (artifacts in the 5.8 /im image 
prevented using this technique at that wavelength). The PSFs we used were obtained from 
the observation of bright stars (Vega, e Eridani, Fomalhaut and e Indi) as part of the Spitzer 
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Space Telescope GTO program 90, and are available at the Spitz er Science Cente r web sit j^. 
The detailed descripti on of this fitting procedure is available in iMarengo et al\ ( 120061 ) and 



Schuster et all ( l2006bl ) . The photometry is listed in Table [21 



2.3. Ground-Based IR Photometry 



We obtained ground-based broadband 1 to 12 fim photometry in August 2000 at the 
University of Minnesota's (UMN) O'Brien Observatory (OBO) with a single element Bolome- 
ter using a standard chop-nod technique for background sky subtraction (beam diameter ~ 
27^\ lGehrz et a/.lll974j . ll992l : lGehrzlll997l ). The observations and photometry are summarized 
in Tables [1] and [2l The absolute calibration of the instrument was detern iined using a Lyr 
(Vega) and /? Peg as primary calibrator, using the procedure described in iGehrzl (119971 ). 



2.4. HST/WFPC2 PC Photometry 



Multi-wavelength images of NML Cyg were obtained i n September 1999 w ith the WFPC2 
Planetary Camera on-board the Hubble Space Telescope (IBiretta et a/.ll2000l ). These obser- 
vations aimed to investi gate the star's CS ne bulosity are summarized in Table [1] (images 
previously published by ISchuster et a/.l l2006al ). NML Cyg was observed with broadband 
Johnson-Cousins V and R filters as well as narrow band Hq, filter. Prior to co-addition, the 
images were processed with the Space Telescope Science Institute's standard calibration ref- 
erence files. The images were combined with the IRAF/STSDAS software package DITHER, 
which uses the Drizzle algorithm to recover imag e resolution from the pixel response of the 



came ra while preserving photometric accuracy (IKoekemoer et all l2000l : iFruchter fc Hook 



20021 ). Multiple, dithered exposures allowed removal of cosmic rays, bad pixels, and other 
effects during co-addition. The HST Photometry is reported in Table [2J 



3. NML Cyg's Bolometric Luminosity and 
10 fiui Silicate Dust Absorption Feature 

Table [2] summarizes our photometric observations of NML Cyg along with archival IRAS 
Point Source fluxes. These observations are plotted on the \F\ vs. A SED in Figure [21 alo ng 



with reprocessed archival ISO SWS spectra (jjusttanont et a/.lll996l : iKraemer et al. 



2002). 



■^http:/ /ssc. spitzer. caltech.cdu/irac/psf.html 
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NML Cyg's spectrum rises rapidly from optical wavelengths through the IR to peak 
in the 2 to 20 nm ra nge, with a broad far-IR tail. Visible light images presented by 
Schuster et all (l2006al ) show that NML Cyg is almost completely obscured by the CS en- 
velope enshrouding the star, thus confirming significant circumstellar extinction. This is 
evident from the broad 10 /im silicate dust absorption feature. Due to the optically thick 
envelope, this feature is a striking characteristic of the NML Cyg s pectrur n, setting this star 
apart from other cool hypergiants such as fj, Cep (IGehrz et a/.lll970l). S Per (lHumphrevslll974 



apart irom otner cooi nypergiants sucn as u uep llLrenrz et aL.\\i.\Jiw. d rer i.M.unipnrevsllib)Y4j . 
and unpubhshed data), VY CMa JSmith et alkoQlh. VX Sgr (lHumphrev] | l97^, HR 5171 A 



flHumphrevs et a/.lll97lr) . IRC +10420 flJones et a/.lll993[ ) and M33 Var A flHumphrevs et al 



20061 ). where the 10 fim silicate feature is seen in emission from optically thin CS nebula 
shells. 



NML Cyg is also a semi-regular variable with period ~ 940 days (IMonnier et a/.l 119971 . 
and references therein) . The vertical sprea d between data sets i n Figure [2] is well within 
the ~ 50% amplitude vari ation reported bv iMonnier et all (119971 ) from 8 to 13 /xm. Even 
with this large variability, IMonnier et a/.r s results show the 10 /im silicate feature's shape 
is roughly constant over a span of nearly 19 years. In addition, the vertical offsets between 
observations at other wavelengths in Figure [2] are at least in part due to variability. However, 
the degree to which NML Cyg's spectrum varies, in either amplitude or shape, at other 
wavelengths has not been well established. Interstellar (IS) extinction to NML Cyg is also 
high, with an estimated range of Ay ~ 3.7 — 4.6 magnitudes (jLedll970l : [Gregory fc Seaquist 



We calculate NML Cyg's minimum bolometric luminosity by integrating^ the SED in 
Figu reBl obtaining Uni = 1.04± 0.05 x 10^ -(d/kpcY Lq, or 3.15±0.74 x 10^ L© at 1.74±0.2 
kpc (IMassey &: Thompson! 1 1 99 ll . Cyg 0B2 distance modulus 11.2). This luminosity does not 
include a correction for IS extinction, which would add at most a few percent to the total 
since most of the light is in the mid-IR where extinction is low. The uncertainty in this 
estimate does not consider the changes in the source brightness due to its variability. This 
luminosity va l ue is i n good agreement with the estimate of Lboi = 1-13 x 10^-((i/kpc)^ Lq by 



Blocker et al\ feOOll ). Thus, NML Cyg's intrinsic luminosity, a few x 10^ Lq {Mm ^ -9.0) 



is similar to that of other extremely luminous M-type hypergiants. 



Previously, iMorris fc Jural (119831 ) reported a luminosity of ~ 5 x 10^ Lq for NML Cyg 
at a distance of 2.0 kpc. This luminosity was equal to the most lumi nous known M su- 
pergiants in the Milky Way and LMC, Mm — -9.5, ~ 5 x lO^L© (IHumphrevsl Il978 



Humphreys fc Davidsonlll979l ). Our value is 37% lower, with 24% due to the downwardly 



*Using a Rayleigh-Jeans F\ (x X "'fit extension at long wavelengths. 
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revised distance and 13% from our integration of the SED from 0.5 to 100 fim (including new 
photometry). With this revision, NML Cyg appears to be less luminous than the hypergiant 
VY CMa (M4-5e la, 4.3 x 10^ L^;^ at 1.5 kpc, also without correcting for IS extinction), which 
has a comparable SED (jSmith et a/.ll200ll : [Humphreys et a/.l 120071 ). The uncertainty in our 
estimate, and the variability of the source makes it difficult to determine unambiguously 
which of these stars is more luminous. Even with this reduction in luminosity, however, 
NML Cyg remains the most luminous star with spectral type M6 or later known in the 
Milky Way. 



4. Mid— IR Image Analysis 



Figure [T] shows that NML Cyg's CS envelope is clearly extended at mid-IR wavelengths, 
compared to the point-source (PSF) images of 7 Dra, as the diffraction pattern appears 
'filled- in'. As further evidence for the extension and degree of asymmetry. Figure [T] also 
shows the NML Cyg images after over-subtracting the PSF (epoch 1, scaled to give zero 
intensity fiux at the center pixel in the subtracted images). The asymmetric CS envelope 
residuals are apparent, having a NW/SE elongation axis, and position angle in the range 
~ 120° — 150°. Also visible in the PSF-subtracted images is a lopsided tail extending to the 
NW (across the first diffraction ring). The MMT is an Elevation/ Azimuth mount telescope, 
so the background sky predictably rotates in the image plane while tracking an object. Thus, 
it was possible to directly confirm the asymmetric extension by observing NML Cyg twice 
in the same night, at two sky rotation angles separated by ~122° (at 9.8 /im only). The 
extended components rotated by the same number of degrees (not shown here), eliminating 
the possibility that the PSF or other instrumental effects are the cause of the observed 
extension. 



NML Cyg's mid- 



H2O (iRichards et al 



R envelope orientation is similar to the 130° — 150° position angle s 



1974 



Benson fc Mutell Il979l : [Diamond et al\ Il984f ) 



obser ved in the OH (iMasheder et al. 

1993)7 (j ClaussenI |2004| . and private communica- 



tion) masers. [Richards et al\ (119961 ) have suggested that the NW/SE spatial distribution 
of the H2O vapor ma sers may indicate a b ipolar outflow, supported by recent SiO ground 
state observations by iBoboltz fc ClaussenI . It is also possible that the maser emission is 
tracing an asymmetric, episodic outflow that may be reminisce nt of the arcs and other 
struc t ures seen in the circurnstellar nebula surroundi ng VY CMa (jSmith et al\ I2OOII : ISmith 
20041 : iHumphrevs et a/.ll2005l . I2OO7I : I Jones et a/.ll2007h . 



Earlier observations by iHinzl (120011 ) also show NML Cyg's extended envelope using 
MIRAC3 with the BLINC module in nulling interferometer mode on the MMT. These ob- 
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servations were made prior to the installation of the adaptive optics system. iHinzl 's observa- 
tions indicated a CS envelope extended beyond 0."28 at 10.3 /im, with nearly equal flux in 
the extended envelope relative to the unresolved structure (see his Eq. 7.5 and Table 7.2, p. 
89-95). The better resolution, and Fourier u, v coverage, of the AO images in Figured] reveal 
at least two structures in the CS nebula, the main envelope aligned with maser observations 
and an asymmetric component aligned more closely with the direction towards Cyg 0B2. 



4.1. Parameterizing the CS Envelope Image Intensity 

To take full advantage of the exceptional stability and contrast in the MMT/AO MIRAC3 
high-angular resolution images, we used deconvolution to more clearly reveal details seen in 
the Figure [1] images of NML Cyg's CS envelope. Figure [3] shows o ur NML Cyg images af 



ter d e convolution w ith the PSF calibrator (using IRAF task LUCY, iRichardsoru 1 19 72l : iLucy 



1974 : ISnyderl Il99d ). The deconvolution reveals the complex CS envelope structure, high- 
lighting the elongated envelope extended towards the NW direction. The two contours in 
the Figure are set at the source 1/2 and 1/10 maximum flux level. The inner contour shows 
that the core of the emission is elongated along a SE - NW axis, while the external part 
of the circumstellar emission is lopsided towards the NW quadrant. The main axis of the 
core emission has a PA of ~ 140° (counter-clockwise from north) in the 9.8 and 11.7 /im 
deconvolved images. This angle is similar to the alignment of the NML Cyg OH, H2O and 
SiO masers, also at ~ 130° - 150°. The PA is slightly smaller (~ 120°) in the 8.8 /zm image, 
possibly due to an image artifact in direction of the MIRAC3/BLINC chop. 

To better characterize the shape and spatial scales of the NML Cyg circumstellar en- 
velope, we flt the images in Figure [3] with two 2-D Gaussians, one fltting the core emission 
(which we call the "core envelope" Gaussian) and one fltting the NW extended low level 
structures (the "outer envelope" Gaussian). Note that this fltting process is aimed only to 
flnd a practical way of describing the changing colors of the circumstellar envelope at differ- 
ent distances from the center. The fltting process or results do not imply that the individual 
Gaussians represent individual physical structures in the circumstellar envelope. The core 
Gaussian has a width and orientation analogous to the width and position angle of the 1/2 
height contour, and is centered on the star. The second Gaussian flts the outer envelope 
emission, and is offset in the NW direction as shown by the 1/10 contour level asymmetry. 
The best flt parameters were determined by minimizing the residuals after subtraction of 
the two Gaussians from the deconvolved images. We start by fltting the flrst Gaussian to 
the image core, then the second Gaussian is flt to the core-subtracted image. The process 
was then iterated to minimize the residuals after subtraction of both Gaussians. Table [3] 
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lists the best fit parameters for eacli wavelengtli, including the size parameter a and relative 
flux (in % of the total flux) of each component, the aspect ratio a/h between the major and 
minor semi-axis and the position angle PA of the core component and the offset of the outer 
envelope with respect to the core. Figure H] shows the radial profiles of the deconvolved 
image and best fit Gaussians along the cut at PA = 140°. Note how the profiles show the 
large asymmetry between the SE and NW directions. While the SE profile is fit well by the 
two Gaussians, some residual flux above the fit is still present in the NW direction. 

The agreement in position angle and size between the core envelope emission and the 
NML Cyg maser observations suggest that this emission arises fro m the NML Cyg's win d 



at, or just interior to, the photo-dissociation front(s) described in ISchuster et al\ (l2006al ). 
The relative flux contribution from the outer envelope component increases from 8.8 to 9.8 
to 11.7 jivci as we are resolving the optically thin warm dust {Tdust ~ 200 - 500 K) farther 
out in the envelope. 



4.2. The Asymmetric NW Excess 

To test the accuracy of the deconvolution parametrization in image space, we have 
subtracted the Gaussian fit, convolved with the PSF reference (7 Dra), from the NML Cyg 
mid-IR images at each wavelengths. The residuals are shown in Figure O The grid-like 
pattern is due to the spatial sampling and linear interpolation of the image on the detector 
grid, and is the limiting factor in the accuracy of the fit. This pattern may be reduced 
by increasing the number of unique offsets in the dithered observations and by applying 
higher order interpolation in the reconstruction of the image on a finer pixel grid. The 
uniform distribution of this residual pattern and the relatively low flux level in the center 
is indicative of the good quahty of the deconvolutions and fits, and also the stability of the 
MMT/AO PSF. 

The figure also shows a clear excess residual in the NW quadrant. The excess peak 
brightness occurs around ~ 0."3 to 0."5 (~ 520 — 870 AU) from the central star, with a 
broad tail extending toward Cyg 0B2. This excess cannot be reduced by adjusting the 
position and shape of the outer envelope Gaussian, constrained by the fit in the other 3 
quadrants. The shape and or ientation of thi s excess are reminiscent to the lopsided emission 



seen in the HST images from ISchuster et al\ (j2006al . their Figures 2 and 6). The asymmetric 



excess is separated from the star in Figure [5] at the spatial resolution limit, i.e. it is resolved 
according to the Rayleigh crite rion. It lies just out side the asymmetric dusty reflection 



nebula seen in the i/^ST images (ISchuster et a/.ll2006al ). We suggest that the excess emission 



is likely from warm dust that is externally heated/destroyed by UV radiation from Cyg 0B2. 
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In the next section, this hypothesis will be tested by directly measuring the flux from this 
excess. 



5. Discussion: CS Dust Heating and Destruction 

In Table H] we list the total flux of the source, the fraction of the total flux in the core 
and outer envelope, the total flux of all the fit residuals shown in Figure [5|, and the excess 
residual flux in the NW quadrant (residual flux in the NW quadrant, minus the residual 
flux in the opposite SE quadrant). In Figure [6] we plot the [8.8]-[9.8] and [9.8]-[11.7] colors 
for each CS nebula component. Both the core envelope colors are bluer than the total for 
the source. This core component likely represents the collective flux from where the silicate 
absorption is deepest, and thus characterizes the geometry of the dust responsible for the 
optically thick feature seen in the star's spectrum. NML Cyg's stellar radius is probably 
of the order of ~10 AU, typical for the largest red supergiants, and therefore the material 
comprising the inner CS envelope {a ~ 0."12, 200 AU) has a size scale of approximately 20 
stellar radii. The outer CS envelope colors are redder than those of the core envelope and 
the total source. This indicates dust with decreasing temperature and density. 

In contrast, the colors for the asymmetric excess are significantly different, being blue in 
[9.8]-[11.7] and red in [8.8]-[9.8]. These colors can be explained by the presence of optically 
thin hot dust with the 9.8 fim silicate feature in emission. The fact that the asymmetric 
dust emission excess appears only on the side facing Cyg 0B2 is consistent with an increase 
in temperature resulting from the external UV radiation. The density is also likely reduced 
due to dust destruction. Since the asymmetric excess is shielded from the central star by 
the optically thick inner CS envelope, Cyg 0B2's external energy is required to account for 
changes in the grain temperature, size and/or spatial distribution ultimately responsible for 
the silicate feature's inversion. 



In ISchuster et all (l2006a] ) we showed that the UV flux from the Cyg 0B2 association 
is strong enough to significantly heat the dust grains in the NML Cyg envelope, possibly 
destroying the smallest grains. The excess mid-IR flux we observe in the NW quadrant may 
be the direct effect of this heating and destruction. Its presence strongly suggests that the 
circumstellar envelope of NML Cyg is indeed shaped by the interaction of its wind with the 
radiation field coming from the Cyg 0B2 association. The direct detection of this physical 
interaction between NML Cyg and the Cyg 0B2 association is further confirmation that the 
source is indeed in th e neighborhood of Cyg 0B2. This supports the distance determination 



of 1.74 lb 0.2 kpc by iMassey fc Thompson! (jl99ll ) for this star as one of the most reliable 



distances for an evolved massive star near the upper luminosity boundary. 
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6. NML Cyg's Complex Circumstellar Environment 



Figure [7] shows tlie NW excess in our 9.8 um deconvolve d image compared to the scat- 
tered light HST/WFPC2 F555W image from ISchuster et all ( l2006al ). As discussed in Sec- 
tion 5, our interpretation for this mid-lR excess is another strong indicator that the UV radi- 
ation from the hundreds of massive young OB stars within Cyg 0B2 is disrupting NML Cyg's 
post-MS wind through photo-ionization, photo-dissociation and grain destruction. The mod- 
els for these processes explain the inverted HIl region's existence and shape, the asymmetric 
nebula seen with HST, as wel l as the externally he ated dust visible in our mid-IR images. 



The in verse photo-ionization (IMorris fc Jural Il983l ) and photo-dissociation (ISchuster et al. 



2006al ) models for this interaction have assumed a radially symmetric, constant velocity 
hydrogen gas wind, i.e. density p cx r~^. The agreement between theory and observation 
suggests that this steady-state, uniform outflow is a overall a good approximation of NML 
Cyg's stellar wind. 

However, the asymmetric mid-IR image reveals that NML Cyg's dusty CS envelope 
does has a more complex underlying structure. One may speculate that the identified struc- 
tures may represent the integrated signal from many arcs, loops, bipolar outflows and other 
3-dimensional structures below our angular resolution that may be the result of episodic 
an d asymmetric mass- los s , much like VY CMa (see Figure 3, ISmith et al\ l200ll as well 



as 



Humphreys et a/.l 120071 : 1 Jones et a/.l 120071 ). Asymmetric structures arising from episodic 



mass-loss may also modulate the maser and IR emission as they pass through the photo- 
dissociation and grain destruction regions. With even better angular resolution and sensitiv- 
ity (perhaps with J VFST and/or the next generation of ground-based nulling interferometers 
like LBTI and KNI) it should be possible to observe these fainter, compact structures within 
the larger nebulosity, if they exist. 



Previous 11 /im results from lMonnier et al\ ( 119971 ) , iDanchi et al\ ( l200ll ) and [Blocker et al. 



(I2OOII ) have suggested the presence of multiple, concentric density enhancements (shells) 
surrounding the central star superposed on a wind, characterized by an azimuthally 
symmetric, but non-steady-state dust outflow, i.e pdust 9^ f'"^- This, in turn, implied a time 
dependent mass-loss, possibly in episodic/periodic events such as a 'superwind' phase. The 
images presented in Section H] show evidence for a complex distribution in the CS dust, but 
also reveal an asymmetric excess aligned to the NW (towards Cyg 0B2). The differences with 
these earlier results are not surprising given the better resol ution and Fouri er m, v coverage 
in the MMT/AO images. However, it should be noted that iMonnier et al\ acknowledge in 
their conclusions that deviations from spherical symmetry, particularly emission from dense 
or clumped material 350 r nas from the centra l star, would necessarily change the meaning of 
their fits (and likewise for lDanchi et a/.ll200ll ). In fact. Figure [3] and H] show, through direct 
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imaging, the presence of warm dust emitting throughout mid-IR wavelengths precisely at 
this distance from NML Cyg. Moreover, this emission is concentrated at position angles 
ranging from ~ —30° (CS env elope) to ^ —6 0° (d ust facing Cyg 0B2), or equivalently at 



alignments of ~ 120° - 150°. iMonnier et a/.r s and iDanchi et alX s journals of observations 



show that they made eight observations with resolution substantially (at least 20%) better 
than ~0."35, and these observations were made with the interferometer aligned along posi- 
tion angles between 102° and 140° (see each author's Table 1). Thus, their observations just 
happened to be oriented with the NW excess emission resolved in our images. It is possible 
that the most external circumstellar shell inferred from these interferometric observations 
may in fact be better explained by the asymmetric excess in our 9.8 /im image. 



Using a combination of Keck aperture masking and IOTA interferometry, iMonnier et al 

( 120041 ) found at 2.2 /im an elongated structure with a diameter of ~ 50 mas. The orientation 
of this structure is orthogon al to the PA of t he inn er component of our images, represented 
by our "core" Gaussian fits. IMonnier et al\ (120041 ) interpret this structure as an equatorial 
enhancement in the dust envelope, orthogonal to NML Cyg's maser outflow. 

Our mid-IR images, when analyzed in combination with the previous observations dis- 
cussed above, clearly demonstrate the complexity of NML Cyg circumstellar environment. 
Any model of NML Cyg that attempts to explain the observations must thus take into ac- 
count the observed asymmetries as well as the physical interaction between NML Cyg's wind 
and the UV radiation from Cyg 0B2. In particular for NML Cyg, high-angular resolution 
(< 0."3) combined with more complete coverage of the Fourier m, v plane are crucial to 
further investigate asymmetries in its complex CS environment. 

NML Cyg is part of our larger program to study the circumstellar environments of cool 
stars near the empirical upper luminosity boundary in the HR diagram with evidence for 
high-mass loss and observed instabilities. As such, NML Cyg was perhaps the best candidate 
to have an extensive CS nebula, possibly like VY CMa. However, NML Cyg's CS nebulosity 
is more co ncentrated ar o und t he star and appears quite different compared to VY CMa, as 
imaged by ISmith et al\ (120011 ). Even though NML Cyg is heavily obscured, the extent of 
the CS nebula is much less than the almost 10" nebula surrounding VY CMa (> lO'* AU at 
VY CMa's distance of 1.5 kpc). It is possible that the more distant material in NML Cyg's 
circumstellar nebula has been largely dissipated by the winds and radiation pressure inside 
the Cygnus-X super bubble and is below our detection limits. If NML Cyg were not located 
in such close proximity to Cyg 0B2, it might show a much more extende d nebula comparabl e 
to VY CMa. However, VY CMa's mass-loss rate, 2 - 4 x 10"^ Meyr'^ JPanchi et a/.lll994j ). 
is 3 to 6 times higher than for NML Cyg and this difference may also have significant bearing 
on the more extensive circumstellar nebulosity. 
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7. Conclusions 

Our sub-arcsec angular resolution AO mid-lR images of NML Cyg provide a new under- 
standing of the complex geometry and physics of the circumstellar environment of this high 
luminosity cool hypergiant. By spatially resolving the optically thick dusty envelope, we 
directly image the structures responsible for the creation of a deep 10 /im silicate absorption 
feature. This structure, which follows the same orientatio n of NML Cyg' s mase r outflow, is 



orthogonal to a near-IR equatorial enhancement found by iMonnier et all ( 120041 ) . 



By analyzing the mid-IR colors of structures located at increasing distance from the star, 
we observe a trend in which the optical depth of the dust decreases in the outer parts of 
the circumstellar envelope. For the first time we isolate an asymmetric excess, at a distance 
of ~ 520 to 870 AU, NW from the star, with colors consistent with the emission of hot, 
optically thin dust. We interpret this emission as the signature of the interaction of the 
NML Cyg circumstellar envelope with the strong UV fiux generated i n the nearby Cyg 0B2 



association. This interaction was predicted in our previous paper ([Schuster et al\ l2006al ) 
to explain the shape of a the inverted photo-dissociation region we imaged with the HST 
at optical wavelength. Our new mid-IR observations strongly support the validity of our 
previous results and the model we proposed for their explanation. 
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his VLA observations. NASA provided support for this work. Observations reported here 
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Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of 
Technology under a contract with NASA. Support for RDG was provided by NASA through 
grants 1215746 and 1256406 issued by JPL-Caltech to the University of Minnesota. This 
work is based in part on observations made with the NASA/ESA Hubble Space Telescope, 
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Table 1. Observations 



iVliVi 1 IViixxAL^o 


A 


A \ 


T-. 

Airmass 


= 

hjxposures 


July id, 2006 U i 


(fim) 


(fim) 




no. X sec 


NMT, Cvtr 


8 80 


88 


1 10 


13x20 


9^80 


0.98 


1^18 


11x20 




9.80 


0.98 


1.02 


2x20 




11.70 


1.12 


1.14 


11x20 


7 Dra (calibrator) 


8.80 


0.88 


1.08 


15x30 




A QQ 
U.OO 


1 1 Q 

l.io 


io X oU 




9.80 


0.98 


1.07 






9^80 


0^98 


1^23 


36x30 




11.70 


1.12 


1.09 


20x30 




11.70 


1.12 


1.28 


20 X 30 


Spitzer/IRAC 


A 


AA 


AOR 


Exposures 


July 27, 2004 UT 








no. X sec 


NML Cyg 


3.550 


0.750 


6588416 


3x10.4'' 


4.493 


1.015 




3x10.4'^ 




5.731 


1.425 




3xl0.4'» 




7.872 


2.905 




3x10.4'' 


OBO MN Bolometer 


A 


AA 


Filter 


Exposure 


Aug 10, 2000 UT 


(/im) 


(/im) 




NML Cyg 


1.250 


0.200 


J 


n/a 


1.653 


0.297 


H 




2.340 


0.500 


K 






3.647 


1.152 


L 






4.900 


0.700 


M 






10.925 


6.730 


N 






7.908 


0.755 








8.808 


0.871 








9.803 


0.953 








10.273 


1.013 








11.696 


1.110 








12.492 


1.157 






HST/WFPC2 


A 


AA 


Filter 


Exposures 


Sept 16, 1999 UT 


((um) 


(nm) 




no. X sec 


NML Cyg 


0.4293 


0.0473 


F439W (B) 


6x500'' 


0.5337 


0.1228 


F555W (V) 


20, 100, 4x400 




0.6564 


0.0022 


F656N (He) 


20, 2x260 




0.6677 


0.0867 


F675W (R) 


0.5, 10 
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Note. — Exposure time not applicable to UMN Bolometer since the in- 
strument measures changes in Voltage to obtain instrumental magnitudes 
and signal-to-noise. 

''Saturated. 

''No detection. 
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Table 2. NML Cyg Photometry 



Instrument Filter 


A (/im) 


r lux \jy J 


m — 


HST/WFPC2 








Sept 16, 1999 UT 








V 


0.5337 


2.377e-04 


0.003e-04 


Ha 


0.6564 


2.395C-03 


0.014C-03 


R 


0.6677 


6.124e-03 


0.029e-03 


OBO MN Bolometer'' 








Aug 10, 2000 UT 








J 


1.250 


21.5 


1.4 


H 


1.653 


131 


0.7 


K 


2.340 


333 


1.0 


L 


3.647 


1287 


4.2 


M 


4.900 


2121 


10 


N 


10.925 


4325 


58 




7.908 


4522 


90 




8.808 


4124 


60 




9.803 


3881 


43 




10.273 


3905 


71 




11.696 


5020 


120 




12.492 


5366 


84 


Sj,'ifzer/mAC PSF fittiua 








July 27, 2001 I T 








3.6 


3.550 


1150 


160 


4.5 


4.493 


1670 


240 


8.0 


7.872 


4160 


830 


MMT AO MIRACS/BLINC" 






July 23, 2006 UT 










8.80 


3735 


63 




9.80 


3780 


160 




11.70 


5280 


130 


IRAS PSC<= 








1983 










12 


5580 


560 




25 


3990 


400 




60 


1030 


100 




100 


335 


34 
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Note. — Photometry not color corrected. 

''Observational errors exclude typical 5 — 10% systematic uncer- 
tainties. 

''NML Cyg's mag is near the limit for IRAC PSF fitting (see 
ISchuster ei af.ll2006b[) . 5.8 fim PSF fitting was not rehablc. 

"^Fluxes from the IRAS Point Source Catalog rejects - from 
the InfraRed Science Archive: http://irsa.ipac.caltach.edu, data 
tag: ADS/IRSA.Gator#2007/1009/122926_27251. 10% flux un- 
certainty assumed. 



Table 3. CS Envelope Component Parameters 



A (//m) 


Comp. 


ARA (mas) ^ 


ADec (mas)^ ^ 


Flux (%)<= 


a (mas)'^ 


a/b 


P.A. 


8.8 fim 


















core envelope 








74.1 


91.0 


1.09 


120 


9.8 Mm 


outer envelope 


14 


87 


25.9 


171.7 


1.00 




















core envelope 








70.4 


124.0 


1.20 


138 




outer envelope 


44 


105 


29.6 


176.5 


1.00 




11.7 
















core envelope 








65.6 


113.1 


1.02 


139 




outer envelope 


45 


98 


34.4 


171.7 


1.00 





^Position offsets are relative to core component, +ARA is to the West 

''The errors, including fit uncertainty and field rotation during the observation, is ~ 5 mas. 

'^Percentages are relative to total fiux. Uncertainties better than 2%. 

dpWHM = 2v^2Tn2 X a 
"Uncertainty of 10 mas or better 

^ Major axis (a) position angle is measured in degrees counter-clockwise from North 
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Table 4. CS Envelope Component Fluxes 



A (/im) Comp. Flux (Jy)*^ 



8.8 /im 




Total Flux 


3735 ± 63 


Core Envelope 


2540 ± 67 


Outer Envelope 


888 ± 23 


Total Residual 


307 


NW-SE Excess 


95 

49.5 


Vega 


9.8 Aim 




Total Flux 


3780 ± 160 


Core Envelope 


2422 ±114 


Outer Envelope 


1018 ± 48 


Total Residual 


343 


NW-SE Excess 


121 ± 4 


Vega 


40.2 


11.7 /im 




Total Flux 


5280 ± 130 


Core Envelope 


3167 ± 101 


Outer Envelope 


1660 zt 53 


Total Residual 


453 


NW-SE Excess 


170 ± 11 


Vega 


28.4 



Note. — Photometry not color corrected. 

^Observational errors exclude typical 5 — 10% 
systematic uncertainties. 

''Negative error includes uncertainties in sub- 
tracting the vertical chop bleeding. 
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Fig. 1.— 8.8, 9.8 and 11.7 /xm MMT/AO MIRAC3/BLINC images of NML Cyg's dusty CS 
nebula. Top row - Images of NML Cyg, no post-processing. Middle row - PSF reference star 
images (7 Dra, epoch 1). Bottom row - NML Cyg after over-subtracting the PSF images 
above (scaled to peak brightness to give zero flux). NML Cyg's CS envelope is clearly broader 
than a point-source, extended with a NW/SE orientation. There is also an asymmetric excess 
on the NW side. The spot in the lower right corner is a PSF ghost. The images are displayed 
with a square root intensity scale. 
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Silicate Feature 



I m i 

rTTTl 
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Fig. 2.— NML Cyg SEP from photomet ry in Table [2], as well as ISO SWS spectra 
(jjusttanont et al. 1996 : Kraemer et al\ [2OO2I ) . The 10 /xm silicate feature is seen in ab- 
sorption, and is a result of the optically thick CS envelope enshrouding the star. There 
is some normalization offset between the ISO spectra and the other observations. At least 
some of this difference is presumably from the variation in the star's spectrum over its 3 yr 
period. The ground-based observational errors have typical 5 — 10% systematic uncertainties 
(only calibration uncertainties are plotted here). NML Cyg's extreme brightness led to large 
IRAC PSF fitting errors. The star's high luminosity may have led to saturation in the IRAS 
observations, resulting in the offset with the ISO spectra. 
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NML Cyg deconvolved 
11.7 /xm 




Arcsec E of Star 



Arcaec E of Star 
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Fig. 3.— MMT/AO MIRAC/BLINC images of NML Cyg after deconvolution with the 
PSF reference star (7 Dra, epoch 1). The inner contour is set to the source 1/2 maximum 
flux level, and shows that the core of the emission is elongated along the SE - NW axis 
(~ 140°, dashed line). The outer contour is set to 1/10 of the source maximum flux level, 
and shows that the external part of the circumstellar emission is lopsided towards the NW 
quadrant. The images are displayed with a square root intensity scale. The arrow indicates 
the direction of the Cyg 0B2 association. 
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Radius [arcsec] Radius [arcsec] 

Fig. 4. — Solid line — Radial profile of the NML Cyg deconvolved images in the SE (left) 
and NW (right) directions. The profile is derived along the dashed line plotted in Figure [3l 
averaged on a 60° wedge. Dashed line — Radial profile of the 2-Gaussian fit in the same 
directions. Dotted lines — Radial profile of the individual Gaussian best fit components. 
Note how the Gaussian fit describes well the profile in the SE direction, but a significant 
excess is left in the NW direction. The deconvolved image profiles are normalized to the 
peak intensities, and the Gaussian profiles have been rescaled accordingly. 



- 28 - 



NML Cyg residuals 
8.8 /im 
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NML Cyg residuals 
9.8 fim 




NML Cyg residuals__ 
11.7 /xm 




Arcsec E of Star 
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Fig. 5. — Residuals after subtraction of the 2-Gaussian fit, convolved with the PSF reference 
(7 Dra, 1st epoch), from the NML Cyg MMT/AO MIRAC3/BLINC images. The residuals 
show a clear asymmetric excess in the NW quadrant, roughly in the direction of the Cyg 
0B2 association. This as ymmetric exces s also a grees with the lopsidedness seen in the 
HST/WFPC2 images from ISchuster et all (j2006al ). External heating by the UV radiation 
from Cyg 0B2 likely causes this excess emission. The bright spot to the North in the 8.8 
lim image is possibly bleeding from MIRACS's vertical chop. Grid-like patterned residuals 
represent the limiting noise due to the image spatial sampling on the coarse detector array 
grid. 
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Fig. 6. — Mid-IR color-color plot for NML Cyg and its CS envelope components from 
MIRAC3 observations. The envelope's core [8.8]-[9.8] and [9.8]-[11.7] colors arc blue as 
compared to the total for the source, indicative of hotter and optically thick dust. The 
outer CS envelope component is redder in both colors, indicating cooler dust. In contrast, 
the colors for the asymmetric excess are blue in [9.8]-[11.7] and red in [8.8]-[9.8], suggesting 
warm and optically thin (9.8 iJ,m silicate feature in emission) dust. The colors are in the 
Vega magnitude system. The larger uncertainty in the NW excess colors reflects the possible 
contribution from MIRAC3 vertical chop bleeding N of the star at 8.8 //m. 
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Fig. 7. — Composite image of the physical interaction taking place in NML Cyg's cir- 
cumstellar environment. Main image - the CS dust heated by Cyg 0B2 is isolated in this 
deconvolved MMT/AO MIRAC3/BLINC 9.8 /im image after CS envelope model subtrac- 
tion. Inset - this HST/ WFPC2 F555W (V ba nd, 0.55 /im) image shows the scattered stellar 
light (reproduced from ISchuster et a/.l l2006a( ). The symmetry axis is shown in the direc- 
tion of the nearby C yg 0B2 st ellar a ssociation. The brightest 22 GHz H2O maser features 
from iRichards et al\ f Fig ure 1. Il996r). an d the molecular photo-dissociation surface projec- 
tion from ISchuster et all (Figure 4. l2006al ) have been superimposed on the image. The white 
mark and large maser near the center indicate the star's location, and maser sizes are roughly 
proportional to their flux. For display the dissociation surface is placed along the outer edge 
of the envelope as seen in the HST image. The 9.8 fim image is displayed with a linear 
intensity scale. The 0.55 /xm image is displayed with a square root intensity scale. 



